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DESIGN OF MULTIVARIABLE FEEDBACK CONTROL 
SYSTEMS VIA SPECTRAL ASSIGNMENT 


By 

Stanley R. Liberty^, Roland R. Mielke^, and Leonard J, Tung^ 

ABSTRACT 

Applied research in the ai'ea of . pectral assignment in multivariable 
systems is reported. A new frequency domain technique for determining the set 
of all stabilizing controllers for a single feedback loop multivariable 
system is described. It is sho\m that decoupling and tracking are achievable 
using this procedure. The technique is illusti*ated with a simple example. 

INTRODUCTION 

# 

This report summarizes the progress of applied research buj.ng conducted 
by the authors and colleagues under the support of NASA Grant NSG 1650. The 
objective of this work is to investigate the applicability of modern control 
theory to the design of multivariable feedback control systems, particularly 
as applied to aircraft flight-control problems. 

Previous research efforts under NASA Grants NSG 1519 and NSG 1650 have 
included an investigation of eigenvaluo/eigenvector assignment procedures 
(refs, 1-3). Early studies focused on an algorithmic formulation of the 
spectral assignment problem by Srinathkumar (ref. 4) , while later studies 
included a geometric formulation of the same problem by Moore (ref. 5), 

Kimura (ref. 6), and Davison and Wang (7,'ef. 7). Based on these theorie.s, 
design procedures have been developed for achieving desired mode mixing 
(ref. S), reducing eigensystem sensitivity to variations in plant parameters 
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(ref. 9), and reducing the effects of actuator niDise an system performance 
(ref. 10). In addition, these procedures have been successfully applied to 
the design of a jet engine controller (I’ef. 11) and an insensitive servo 
control system (refs, 12, 13). A detailed .survey of eignevalue/eigenvector 
assignment design techniques was recently presented by the authors in a 
21-hour short course entitled "A Geometric Approach to Control System 
Synthesis" at the NASA/Langley Research Center (LaRC) , The course was 
extremely well received and demonsti*ated that eigenvalue/eigenvectoi.* 
assignment techniques are accessible to practicing conti’ol engineers. 

Research efforts during the past seven months have focused on 
frequency domain theories utilicing spectral factorisation techniques of 
the Wiener-Hopf class. The Wiener-Hopf theories have been studied by 
sevei'al investigators including Whitbeck (rets, 14-16), and were developed 
in the most generrd form by Youla (refs. 17, 18). The Youla wor’k deals 
with the design of optimal controllers for multi-input-output plants 
imbedded in a single multivariable feedback loop configuration. Two 
significant contributions were made in this work. First, a procedure was 
given for characterising the set of all controller's of a particular class 
which stabilise the overall system; second, a method for selecting an 
optimal controller from the set of all stabilising controllers was presented. 
The optimization procedure wvs based on the standard IViener-Hopf filtering 
problem. Youla' s stabilisation procedure has been further studied and 
generalized by Desoer (ref. 19). Based on this work, a frequency domain 
design procedure for consti'Ucting a stabilizing controller has been 
developed, This pi'ocedure is explained and illustrated with an academic 
example in the next section. 

DESIGN OF A STABILIZING CONTROLLER 

Consider the single feedback loop multivariable control system shown 
in figure 1. In this figiire, P(s) is a proper rational matrix transfer 
function representing the system plant, while C(s) is a proper rational 
transfer function matrix representing the conta'oller. It is desired to 
first characterize the set of all transfer function matrices C(s), such 
that the closed-loop system is stable. Then, a search will be conducted 









over this set o£ stabilizing controllers to locate the particular conti'oller 
which also achieves desired mode decoupling and/or system tracking. The 
following definitions are essential to this development; 

Definition 1 ; A proper rational matri.x is said to be exponentially 
stable if all poles of each matrix elenent are located in the left half of 
the complex plane. 

Definition 3, : Given a proper rational matrix P(s), the expression 

P(s) = right exponentially stable rational 

fraction description (right ESRFD) of P(s]| if both Nj,(s) and ^'^t'e 

proper and e.xponentially stable. In a similar manner, the expression P(s) 

= Dj^”^(s) Nj^Cs) is said to be a left exponentially stable rvational fraction 
description (left ESRPO) of P(s) if both N^^Cs) and Dj^(s) are proper 
and exponentially stable. 

Definition 5 : A right ESRFD of P(s), P(s) = N^,(s) D^/^(s), is said 

to be coprime if there s.xist matrices U^(s) and V^,(s), both proper and 
exponentially stable, such that Uy(s) + Vy(s) D^,(s1 = 1. Similarly, a 

left ESRFD of P(s), P(s) = (s) N^(s), is said to be coprime if there 

exist matrices Uj^(sD and V^(s), both proper and exponentially stable, 
such that Njj^(s) Ujj^(s) + Dj^(s) Vj^(s) = 1. 

The above definitions make possible the following characterization of 
the set ox all stabilizing controllers for the system of figure 1. A proof 
of this theorem is given in reference 19. 

Theorem 1 : Let P(s) = N^.(s) D^,"^(s) be a right coprime ESRFD of P(s) 

with U^,(s) N^(s) + V^(s) D^{s) = 1, and let P(s) = D (s) Nj_(s) be a left 
coprime ESRFD of P(s) with Nj^(s) U^(s) + D^(s) Vjj^(s) = 1- Then, the 
closed-loop system is stable if C(s) is chosen as 

C(s) = [W(s) Nj^(s) + V^.(s)]"^ [-W(s) Dj^(s3 + U^,(s3] 

where W(s) is any proper exponentially- stable matrix such that lw(s) Njj^(s) 
\Cs)i ^ 0* Moreover, C(s) is a proper rational macrix. The transfer 
function matrix relating output to input, T^Cs), is given by 

Ty(s) = N^,(s) [-W(s) Djj^(s) + U^(s)] 
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and the transfer function relating the error signal to input, T (s) , is 
given by 

TgCs) = 1 + TyCs). 

The freedom available in choosing a stabilizing controller C(s) is 
thus characterized by the freedom available to choose W(s). If the only 
design consideration is to stabilize the overall system, then any proper 
exponentially stable matrix WCs) will yield an appropriate controller. 

If, in addition, output decoupling is desired, then a proper exponentially 
stable W(s) such that T^(s) is diagonal is chosen. Similarly, if it 
is desired to tra.ik step inputs, a proper exponentially stable iV'Cs) for 
lim 1 

which s-»-0 exists in chosen. 

Problems in which decoupling and tracking are desired are simplified 
when PCs) is exponentially stable. In this situation, N^j^Cs) = N^,(s) 

= PCs), Djs) = 1, D^.Cs) = 1, U^,Cs) = 0, Uj^Cs) = 0, V^Cs) = 1, and V^.Cs) = 1. 
Then, 


TyCs) -PCs) WCs), 
TgCs) = 1 + PCs) WCs) 


and 


CCs) = -WCs) [1 + PCs) WCs)]"^ 


An example will now be presented to illustrate this technique, 
Consider a two input, two output plant described by the transfer matrix 


PCs) 


PliCs) P12CS) 

PjlCs) P22CS) 


1 

Cs + 2) Cs^ + 4s + 2) 


Cs + 1) Cs + 3) 


Cs + 1) Cs + 3) 


S 


Since PCs) is exponentially stable, stability of the closed-loop system 
is assured if the controller transfer matrix is chosen as C(s) * -W(s) 

X [1 + p(s) W(s)]"^, I'ihere WCs) is any exponentially stable rational 
matrix having appropriate dimension. With this controller, the output 
transfer function is Ty,(s) = -P(s) W(s) and the error transfer function is 
TgCs) = 1 + PCs) WCs), 

If decoupling is desired, WCs) must be chosen such that Ty.C3) 

= -P(s) WCs) is a diagonal matrix. With WCs) given by 


WCs) = 


Wii(s) 
W21 Cs) 


'VisCs) 

W22CS) 


the closed-loop transfer function is 



PilCs) 

P12CS) 


WiiCs) 

W12CS) 

T Cs) = -PCs) WCs) = 





y 

.P21CS) 

P 22 Cs)_ 


S 2 X Cs) 

w^aCi-) 


PllCs) Wii(s) + P12CS) WjiCs) PiiCs) WiaCs) + P12CS) W22CS) 

P21CS) WnCs) + P22CS) Wi 2 (s) PziCs) Wx 2 (s) + P22CS) W22CS) 


Thus, for decoupling. 


-P12CS) 

W12CS) = pYx (s)" "^22 Cs) 


1 

Cs + 1) Cs + 3) 


W22 (s) 


and 


-P21CS) 1 

W21CS) = "^11 = Cs + 1) Cs + 3)' "'saCs). 

With this constraint on WCs), 


PCs)l 

P22CS) 


WxiCs) 


PCs) 

PiiCs) 


W12CS) 


Ty(s5 = - 


0 


(s +^1) (s^'- 3) 1 


(S + 1 ) (s"+ 3) "^22 (s 3 J 


If, in addition to decoupling, it is desired that the system track step 
inputs, then every element of the error transfer function T (s) must 
contain a factor s in the numerator, For the previous selection of 
WCs), 


1 


Cs + 2) 

(s + 1) (s + 3) 


Wii(s3 



TeCs3 = 


0 


1 + 


Cs + 2) 

Cs + 13 (s + 33 


W22Cs3 


Thus, there axe many choices of WCs3 which will result in the desired 
tracking propei-ty. A particulai’ly simple choice is WuCs3 = h'22Cs3 = - 
Then , 


WCs) = - 


2Cs + 1) Cs + 3) 


2Cs + 1) Cs + 3) 


£ 

2 


and the controller ti'ansfer matrix becomes 
5Cs + 1) Cs 3) 


CCs) = 


2 S 


FI) 


2 s (s + f ) 
5Cs + 13 Cs + 5) 


[ 2 s (s .|) 2 s (s .|) _ 


TI\e svstem transfer functions are then given by 
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NljOJ 





3(s + 2) 

2Cs + 1) (s + 3) 


and 


T,(s) = 


£di 


Cs + i) (s + 3) 


3(s 4- 2 ) , , 

2(3 + 1) Cs + 3’)"'^ 


•kill 


(s + 1) Cs + 3)J 


The frequency domain technique demonstrated above appears to be an 
attractive method for control system design and is worthy of further study. 
Investigation of two areas in particular appear promising. First, when 
the above techniques are applied, the set of all stabilizing controllers 
for a given system configuration is obtained. The freedom available to 
the designer is to choose a particular controller from this set by 
specifying a proper exponentiariy stable, but otherwise arbitrary, matrix 
W(s), which is selected to achieve a desired system characteristic such 
as mode decoupling, reference tracking, or disturbance rejection. However, 
a systematic procedure for generating a HCs) to achieve given design 
specifications has not been thoroughly developed. Second, the stabilization 
techniques have been developed only for the single multivariable feedback 
loop configuration. Preliminary studies, however, have indicated that a 
similar formulation may be possible for other system topologies. In 
addition, several classes of problems have been identified in which the 
use of other topologies would greatly enhance the freedom available to the 
designer. 

In addition to being a possibly attractive design procedure, the 
generality of the above approach holds promise of a tie with the 
spectral assignment techniques previously investigated. Investigation of 
these relationships should significantly enhance the applicability of the 
stabilization techniques and improve the spectral assignment procedures. 
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EIGENVECTOR ASSIGNMENT FOR NOISE SUPPRESSION 
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ABSTRACT 


A new method of selecting a multivariable state feedback controller is 
presented. The resulting controller simultaneously realizes arbitrary dis- 
tinct closed-loop eigenvalues, approximi'tes specified modal behavior, and 
reduces the effevOts on the system of actuator noise. The method character- 
izes the subspace associated with each closed-loop eigenvalue in which the 
corresponding eignevector must lie. The vector that is to be approximated 
by this eigenvector is projected onto the subspace. When all the eigenvectors 
have been obtained in this way, the steady-state mean-squared error of the 
state variable of the closed-loop system is reduced local to this set of 
eigenvectors using a gradient search technique. 


I . INTRODUCTION 


Consider a linear, time-invariant, multivariable system described by 
X * Ax + Bu + C, where x, ^ E R^, u E R'”, and the pair (A, B) is controllable. 
In this paper consideration is given to the eigenvalue/eigenvector assign- 
ment problem of selecting a stationary feedback matrix F in the control law 

A 

u ■ Fx such that the eigenvalues of the closed-loop system matrix A = A + BF 
are arbitrarily placed. Throughout this paper the closed -loop eigenvalues 
are assumed to be distinct. Once the eigenvalues are assigned, freedom re- 

A 

maining in selecting the eigenvectors of A is utilized to approximate 
specified modal brhavior and to reduce the effects on the closed-loop sys- 
tem of the actuator noise, 

A 

It has been shown by Moore [1] that Che eigenvalues of A can be arbi- 

A 

trarily assigned, but that each eigenvector of A must lie in a subspace of 
r” determined by the corresponding eigenvalue. However, despite the re- 
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stvlctions imposed on the eigenvectors, there remains some freedom in select- 
ing the eigenvectors. This freedom can be used to affect mode mixing [2, 3], 
reduce system sensitivity to plant parameter variations [4], arid as is 
shown in this paper, reduce the effects of actuator noise. 

An initial eigenvector assignment is made by successively projecting 
each vector from a specified set onto the corresponding subspace. This 
yields a set of eigenvectors each of which is the best approximation in a 
least-squares sense to the corresponding specified vector. The eigenvectors 
obtained are c:cpressed as linear combinations of basis vectors for the sub- 
spaces. The steady-state mean-squared error of the state variable is ex- 
pressed as a function of the scalar coefficients in these linear combinations, 
and is reduced local to the initially assigned eigenvectors by a gradient 
search procedure over the scalar coefficients. 

II. INITIAL EIGENSYSTEM ASSIGNMENT 

Let ..., be a self-conjugate set of distinct complex numbers. 

Then it is known that there exists a real-valued matrix F such that 

(A + BF) v^ =» ^i''i» ^ 

if and only if the following three conditions are met for i » 1, ..., n: 

(i ) Vectors v^ are linearly independent in c'^; 

(ii) V, * V.* whenever X. ■* X.* ; 

N / i j ^ j > 

(iii) v.£ span CN, }, where the columns of K, * 

1 

ker[X^I-A:B] . Here, is partitioned compatibly with [X^I-AiB], 

If {pj^ p^} is a set of vectors, and it is desired that the ith eigen- 

vector approximate p^ , i = 1, . . . , n, then the closest approximation in a 
least-squares sense is obtained by projecting p. onto span {N, }. This 

1 A , 


N, 


form a basis for 
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projection is expressed as 


^i 


a ) 


where b. Is the ]th column oj Nj and .... ) is e set ot scalars, 




or equivalently, 


'-1 ■ 


( 2 ) 


m 

where . . . ct^ ) , The p^ vectors are projected succsssively 




and when all have been projected, F is calculated frotn 

•1 

t IV 

. n'* 
n 




(3) 


Where V - [v. 


V ]. 

n'' 


III. NOISE REDUCTION 


Consider a system with actuator noise described by 

X « Ax + Bu + 5 , (A) 

where C is assumed to be aero mean white noise characterized by Che 
following equations 

E{5(c)> - 0 , (5) 

E{^t)C'^(t)} - E 6(t-T), (6) 

E{x(o)> - Xq , (7) 

E{[x(o) - x^][x(o) “ x^]'^) - P^, (3) 

and eC[x(o) - x^jc'^Ct)} - E{5(t)[x(o) «• - 0. (9) 
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The effect of actuator noise on the response of the system can be reduced 
by proper choice of the feedback gain matrix F. Noise reduction is made 
in the sense that the steady-statr "an-squared error of the state vector 
x(t) is reduced; that is, 

lim E{ 1 |x(t) - x(t) I 1^} 

t*t<» 

is reduced where x(t) ■ E{x(t)}. 

To see how this is done, it is helpful to note that the mean-squared 
error is equal to the trace of the state covariance matrix, W(t). Thus, 

lim E{||x(t) - x(t)l|^} - tr W (10) 

t-w 

where W ■ lim W(t). Moreover, if *(t) is the state transition matrix 
ss 


for the closed-loop system, then 

W(t) - 4>(t-T)5<J>^(t-T)d^ + $(t)Po$'^(t). (11) 


Diffexentiating (11) with respect to t yields 

W(t) - A W(t) + W(t) + H. (12) 

A 

If A has eigenvalues in the left-hand half -plane, then lim W(t) exists 

t-HO 

• 

and therefore lim W(t) ■ 0. Thus, equation (12) reduces to the Lyqpunov 

t4« 

stability equation 

^ ^rp 

AW +WA +H-0 (13) 

ss ss 

whose solution is obtained from 


vec W 

ss 


~[I®A + A®1] ^vec 


(14) 


where 


vec [a, • • ’ a 1 

'-1 n '' 
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and 


A® B 


“ll® *12 

*21® *22®‘"*2n® 


a , B a , B ••• a B 
ml m2 mn 


for A ■ [a. . 

Ij mxn 

Since is a function of it is a function of the eigenvectors 
and, therefore, of the scalars in (1). Thus, 

lim E{jjx(t) "■ x(t)ll^) ■ tr[W ] can be reduced local to the initial eigen- 
vector assignment through a gradient search over these scalars. The gra- 
dient, ■ V(trW ) , is a vector with components of the form 

5 "[“aal 


■ tr 


3 X 


ij 


vec ^ |a + B^jV“^Ol)a-^ vecnj (15) 

mm» mm 


where a" + and B^^ is a matrix whose only nonzero column is the 

til ^ •.J, 

i^ column, which is - S) . The symbol "vec " indicates the oper- 

2 

ation of "unstacking" the n -component vector into an n ^ n matrix, The 
steady-state mean-squared error is reduced interatively until an acceptable 
trade-off is reached between obtaining eigenvectors that closely approximate 
the desired vectors and obtaining a small steady-state error. 
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